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Abstract. We show that a certain class of flare models 
for variability from accretion disk coronae are subject to 
beam-plasma instabilities. These instabilities can prevent 
significant direct acceleration and greatly reduce the vari- 
able X-ray emission argued to arise via inverse Compton 
scattering involving relativistic electrons in beams and soft 
photons from the disk. 
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1. Introduction 

The origin of fluctuations in the emission from Active 
Galactic Nuclei (AGN) and binary X-ray sources is an im- 
portant and long-standing problem. One frequently con- 
sidered possibility employs flares in the coronae around 
accretion disks to produce rapid energy release, particle 
acceleration and radiation (e.g., Galeev, Rosner & Vaiana 



197S ; Kuperus & Ionson 1985; for a review, see Kuijpers 
1995 ). These models usually build upon our understand- 



ing of solar flare physics. 

A particular model of this type has been proposed by 
de Vries & Kuijpers (1992; hereafter dVK), and was specif- 
ically applied by them to X-ray variability of AGNs. Their 
model is an elaboration on typical flare scenarios, in that, 
as usual, the source of energy is stored in magnetic fields 
in coronae of accretion disks. They estimate the power re- 
leased in flares in a radiation pressure dominated corona, 
which they stress is a different environment from the gas 
pressure dominated solar corona. They argue this leads to 
a situation where beams of relativistic electrons are pro- 
duced in the corona and then lose essentially all of their 
energy through inverse Compton scattering on UV disk 
photons before they can stream back to the disk. They 
further argue that these inverse Compton (IC) photons 
produce the X-ray variability seen in Seyfert galaxies, and 
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are able to calculate spectral power-densities in reasonable 
agreement with observations. 

However, the dVK model does not take into account 
other mechanisms that might vitiate some of their key 
assumptions. We note that dVK briefly argue that, par- 
ticularly if radiation pressure dominates the energy den- 
sity in the corona, as is indeed likely around standard 
thin accretion disks (e.g., Shakura & Sunyaev 1973 ) which 
they assume, energy losses through scattering on plasma 
waves are unimportant; then the dominant losses will be 
to IC scattering. However, it is well known that an electron 
beam-plasma system is often susceptible to the excita- 
tion of beam-plasma instabilities which usually have large 
growth rates (Sturrock 1964). Here we argue that when 
these beam-plasma instabilities (BPIs) are taken into ac- 
count, the rate of loss of energy by the electrons for the 
accretion disk coronae conditions suggested by dVK is typ- 
ically much higher than the rate of gain of energy through 
direct acceleration by the electric fields, which are pre- 
sumed to arise in reconnection events. Therefore beams 
of electrons usually will not reach the high Lorentz fac- 
tors needed to produce most X-rays by the IC process. In 
many accretion disk models X-rays are usually produced 
through IC scattering of soft photons on h ot the rmal elec- 
trons (e^g. Shapiro, Lightman & Eardley 1976 ; Liang & 
Price 1977 ) . In such a situation beam-plasma instabilities 
are not excited, and only thermal spontaneously excited 
plasma waves should exist. These will have energy densi- 
ties less than the thermal energy density of the plasma, 
which in turn is much less than the radiation energy den- 
sity. In this case, the argument of dVK would be valid, but, 
once they assume a beam is present, then beam-plasma in- 
stability effects must be included. 

2. Growth of Beam-Plasma Instabilities 

The key assumptions of the dVK model arc that: 1) relax- 
ation of magnetic structures efficiently produce relativistic 
electron beams; 2) the particle beam is a mono-energetic 
stream of electrons with an initial Lorentz factor 70; 3) 
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the ambient radiation is from a quasi-infinite disk and can 
be considered as uniform and isotropic, with a radiation 
density u ra( j; 4) the beam is optically thin, so multiple 
scattering of photons can be ignored. Although (3) is an 
approximation, it is a reasonable one, and (4) is certainly 
plausible under many circumstances. But the core of their 
argument hinges on the ability of the neutral sheet in the 
reconnection process to quickly accelerate electrons via a 
direct electric field. During this acceleration process dVK 
claim the equation for the acceleration of an single elec- 
tron suffering 1C losses is 
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where, xi = eE/m e c and \i — ^TU-rad/^n^eC, with all 
symbols having their usual meanings. In that the first 
(positive) term starts out substantially greater in mag- 
nitude than the second (negative) one, acceleration will 
ensue until a limiting Lorentz factor is reached when the 
two terms balance: 
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The electric field is reasonably taken by dVK to be the 
Dreicer value, which we take as: Er> = 67rn p e 3 lnA/(fcsT e ), 
where n p is the electron density of the ambient plasma, 
InA « 20 is the Coulomb logarithm, and all other sym- 
bols have their usual meanings. With typical AGN values 
(np ps I0 10 cm" 3 , T e PS 10 6 K, and T rad ps I0 5 K) they find 
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30. They then conclude that the elec- 



trons will all reach this terminal Lorentz factor before the 
acceleration terminates and the electrons then lose their 
energy against the disk photons providing the background 
radiation field. 

We now show that since a BPI is excited, it will dom- 
inate the energy losses for the beam and actually prevent 
the electrons from reaching the high Lorentz factors cal- 
culated by dVK. Under these circumstances there will be 
very little IC radiation, so that, while a great deal of en- 
ergy may be released through magnetic reconnection, the 
bulk of the energy will pro bably provide heating to the 
corona (e.g., Liang & Price 1977 ) but is unlikely to yield 
the bulk of the X-rays directly through IC emission. 

The dominant growth rate of the BPI depends on the 
relative magnitudes of the bulk velocity of the beam, Vb, 
and the mean thermal velocity in the beam, VTb', under 
some conditions, vt s , the mean thermal velocity of the 
ambient electrons, also must be taken into account. The 
standard formula for the BPI growth rate, valid for yb > 
(rip/nbY^VTb, is our Case 1 (e.g., Mikhailovskii 1974 ) 
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where rib is the beam density, n p is the ambient plasma 
density (here, in the disk corona), and oj pe = 5.47x 10 4 n 
is the plasma frequency in terms of the ambient electron 
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number density in cgs units. The frequency at which this 
mode grows is uj pe (1 - 0.4(ni/np)V3). 

If the beam starts out very slowly, with i% < 
{rip/nbY^VTb, then the "weak" version of the BPI is rel- 
evant, and this is our Case 2 (e.g., Benz 1993) 
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and the frequency at which this dominant mode is excited 
is u) pe . Under the limited circumstances that vt s > i>6 > 
vxb, the "hot-electron" Case 3 yields (e.g., Mikhailovskii 
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where this dominant mode is at a frequency of 
(v b /v Te )uj pe . 

The AGN corona values of dVK for n p = 10 10 cm -3 
and T = 10 6 K, which we also believe are reasonable, will 
be adopted here. There are, however, additional parame- 
ters that must be considered now (basically in lieu of the 
radiation temperature, or u ra d, needed by dVK). First, 
Q = ri b/n p ; for solar flares this value is ~ 10~ 6 - 10~ 4 
(Benz 1993); however, we will bear in mind the possibil- 
ity that this ratio may be higher in this type of radiation 
dominated plasma. We also need initial values of Vb and 
vxb, to determine which of the three Cases defined above 
should be considered. For us to say that a beam actually 
exists we must always demand that Vb > VTb- 

Note that the BPI directly gives the rate of growth of 
an electric field in the plasma, and the energy loss goes as 
the square of the field strength. Then we find that when 
the relativistic effects that arise if the Lorentz factors re- 
ally could become large are included, the rate of change 
of energy of electrons in the beam is, 
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where a = W/E — W/^nbmc 2 , is the ratio between the 
wave energy density, W , and the electron beam energy 
density, E. In order to determine W, knowledge of the 
saturation mechanisms of the wave field are needed. Of- 
ten, in order to avoid a detailed discussion of the satu- 
ration mechanisms, which tend to operate in multiplicity 
in a plasma, the condition of equipartition of energy be- 
tween the waves and t he beam particles is used (Treumann 
& Baumjohann 1997). In that case a is approximated to 
unity, and we consider this situation first. Case 4, where 
a <C 1, and the saturation occurs earlier by trapping, will 
then be addressed. 

In Eqn. (6) we have ignored the IC term appearing 
in Eqn. (1), having replaced it with a generic form of the 
BPI growth rate; the fact that the BPI term is much bigger 
than the xi term for all reasonable circumstances will soon 
become evident. The dominant dependence of upon 7 
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for the first three cases a rises through the replacement: 
n b — > n b / 7 3 (e.g., Walsh |l980| ; Krishan 1999), which ef- 
fectively modifies £, which is defined as the density ratio 
at non-relativistic relative velocities. In Cases (2) and (3) 
we must also write v^/c = (7 2 — l) 1 / 2 /7. 

Now, for Cases 1, 2 and 3, respectively, we have: 
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with Ai — OJ^^utpe ~8x ^ 7 C^i^Jio^ where the com- 
mon notation, X n = X/10 n , has been employed so that 
the physical parameters will be of order unity; 
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with Ai = 0.5(c/z;Tb) 2 C a -'pe 
we have now defined rjb = VTb/c 
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where 



7 (7 3 

with A 3 = C 1/2 (v T e/c)uj pe : 
now, r/ e = vxe/c ~ 0.01. 

Under any of these situations we have %i = 5-2 X 
lO^^ne^ioT^g 1 with our definition of Ed (which is slightly 
larger than that of dVK, thereby only strengthening our 
argument). For any plausible initial value of 7 ~ 1 the 
different dependences of Eqns. (7-9) upon 7 are not im- 
portant. What is important is that Ai,A 2 ,A 3 ^> \x ^ X2] 
i.e., the energy loss term arising from any form of the BPI 
completely dominates over the energy gain term from di- 
rect electric field acceleration. 

We now consider Case 4, where equipartition is not 
established. Under these circumstances, the growth of the 
Langmuir waves for the fastest initial beam situation, Case 
1, is arrested by the trapping of the beam electrons. In this 
case, the ratio a is eventually given by the saturated value 
(Melrose [1986] Krishan |l99|), a = 9/2[n fc /(2n p7 3 )] 2 / 3 , 
and it can be a rather small number that reduces the loss 
rate significantly. This gives a chance for the situation en- 
visioned by dVK to occur. In addition, a initially can start 
out below the saturation value as it arises from thermal 
fluctuations, and thus it could allow an initial thermal run- 
away. The detailed spatial and temporal structure of the 
reconnection sites will determine if this initial acceleration 
can play a significant role. 

In spite of these uncertainties, we can obtain a rea- 
sonable estimate of the influence of BPI in the situation 
where equipartition is not established. We again consider 
all three cases discussed above, but we now include elec- 
tron trapping and assume a to take the saturation value. 
Here the competition between the IC losses represented 
by X2 and the BPI losses represented by the Melrose a 
has to be considered carefully. 

The inverse Compton term increases with 7 whereas 
the a factor modifying the BPI term decreases with 7. 



Thus demanding that the BPI term is smaller than the IC 
term fixes the minimum value of 7 necessary to validate 
the dVK proposal. A detailed calculation yields the results 
for the three cases as follows: case (la), 7™,™ = 54; case 
(2a), 7 m i n = 21; case (3a), 7 mm = 9.16. Thus it is clear 
that only in case (3a), is it likely that the IC term domi- 
nates and hence the dVK proposal is valid. This requires 
rather special conditions for the flare models to work. 

This type of runaway acceleration has been observed in 
the laboratory under specific circumstances which lead to 
a very weak beam plasma instability. In laboratory exper- 
iments, the runaway electrons are observed detached from 
the main body of the plasma, as for example in a stellara- 
tor. If the runaway electrons hit the tungsten aperture, 
they generate X-rays which can be detected. Provided the 
conditions are right, the runaway electrons undergo insta- 
bilities producing plasma oscillations which then couple 
to the ions. This principle is applied in the design of some 
electron tube oscillators (Rose & Clark 1961). Thus the 
runaway electrons can stably propagate under certain cir- 
cumstances, but will be affected by a BPI if they do satisfy 
the conditions for it. These conditions are essentially on 
the velocity of the beam and its thermal spread, as we 
have already discussed for the first three cases above. 

Often it is found that a regime of strong Langmuir 
waves is quickly reached and these waves are further sub- 
jected to modulational instabilities. Thus, different satu- 
ration mechanisms operate at different stages of the de- 
velopment of the instability, depending on beam plasma 
parameters. However, under the circumstances and pa- 
rameters proposed by dVK, the damping is severe. 

3. Discussion and Conclusions 

We thus conclude that the mechanism proposed by dVK 
should not generally work unless much greater densities 
are possible in the coronae at the same time that the tem- 
peratures are lower, since xi rises faster with n p than 
does any form of Tf, p , and declines faster with temper- 
ature. While denser coronae should be available around 
the accretion disks in X-ray binaries, the ambient temper- 
atures will also be a good deal higher, so we cannot suggest 
a physically interesting situation where the BPIs do not 
dominate. If one could somehow begin with very large 7 
values, then the growth rate of the beam-plasma insta- 
bilities are reduced. For Case 1 this does not help, and 
no solutions for large 7 are possible; however, for Cases 
2 and 3, the relativistic decreases in the BPI rates are 
so substantial that high asymptotic 7 values are allowed. 
This is also true for Case 4, where the saturation reduces 
the effectiveness of the BPI; however, even then the BPI 
can prevent much acceleration unless the beam already 
starts with a substantial value of 7 or has such a low den- 
sity in comparison to the ambient medium that it could 
not carry significant power. Moreover, we see no way to 
achieve these initially high 7 values: that is what the dVK 



4 



Krishan et al.: Beam-plasma instabilities 



mechanism was supposed to accomplish, but now appears 
to be incapable of achieving. 

Filamentation, which could produce denser beam frag- 
ments, could play a role by raising £ locally. If any analogy 
can be drawn with solar flares, then the presence of rapid 
irregularities within the Type 3 radio bursts strongly in- 
dicates that the flux tubes are filamenta ry dur ing the ac- 
celeration phase (e.g., Vlahos & Raoult 1995 ). However, 
this possibility is still insufficient to salvage this mecha- 
nism for AGN coronae, since even with £ ~ 1 the ratios of 
^1,2,3/xi > 1- In Case 4, where saturation is important 
in principle, the large value of £ implies that a ~ 1 too 
(for initial 7 ~ 1) so the loss term still would dominate. 

Nonetheless, even with much of the energy going into 
wave turbulence, as we have argued, significant IC emis- 
sion can be possible. This is because (as pointed out by the 
referee) trapping and other nonlinear effects can roughly 
heat the electrons up to kT e ~ e</>, with cj) the electro- 
static amplitude of the waves. Since the energy gain term 
(the first on the RHS of Eq. [1]) is essentially a constant, 
these 'thermalized/trapped' electrons can attain nearly 
the same energy as in the dVK picture. However this en- 
ergy will not be in the form of a beam, as argued by dVK, 
but rather, will be present in an isotropic distribution. 
Then the IC process still works, and one of the points 
made by dVK, that much of the energy is lost by IC hard 
X-rays instead of 'soft' X-rays from material evaporated 
from the disk, can remain valid, as already noted at the 
end of §1. In order to see if the inverse Compton losses 
actually dominate, detailed computations of these effects 
should be undertaken under various circumstances. 

It is well known that in the case of the solar corona, the 
directly accelerated beams should be thermaliz ed wi thin 
a very short time through BPI (e.g., Sturrock 1964 ). In 
the standard picture, this produces Langmuir waves which 
then manifest themselves as various types of radio bursts 
if non-linear effects or transport from faster to slower 
electrons within the beam could dominate (e.g., Vlahos 
& Raoult 1995 ). However, energetic electrons have been 
observed in satellite measurements in near-earth orbit, 
and the outstanding question of the maintenance of these 
beams through their propagation from the sun to the earth 
has given rise to more complex models involving complex 
profiles of the electron beams (Vlahos & Raoult 1995| ). 
Instead of producing X-ray flares via a primary process 
as proposed by dVK, these secondary processes involving 
energy input to the plasma could contribute to variability 
in the radio band. 
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